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We review various modified gravities considered as gravitational alternative for darlc energy. 
Specifically, we consider the versions of f{R), f{G) or f{R, G) gravity, model with non-linear grav- 
itational coupling or string-inspired model with Gauss-Bonnet-dilaton coupling in the late universe 
where they lead to cosmic speed-up. It is shown that some of such theories may pass the Solar 
System tests. On the same time, it is demonstrated that they have quite rich cosmological struc- 
ture: they may naturally describe the effective (cosmological constant, quintessence or phantom) 
late-time era with a possible transition from decceleration to acceleration thanks to gravitational 
terms which increase with scalar curvature decrease. The possibility to explain the coincidence 
problem as the manifestation of the universe expansion in such models is mentioned. The late 
(phantom or quintessence) universe filled with dark fluid with inhomogeneous equation of state 
(where inhomogeneous terms are originated from the modified gravity) is also described. 

I. INTRODUCTION 

The dark energy problem (for recent review see |3jI3) or, why current universe is expanding with the acceleration, 
is considered to be the one of the most fundamental theoretical problems of XXI century. There are various directions 
aimed to construct the acceptable dark energy model. Specifically, one can mention scalar (quintessence or phantom) 
models, dark fluid with complicated equation of state (EoS), more complicated field theories with fermions, abelian 
or non-abelian vector field, string/M-theory, higher dimensions, etc. Nevertheless, despite the number of attempts 
still there is no any satisfactory explanation of dark energy origin. This is understandable having in mind that even 
current values of cosmological parameters are not yet defined with precise accuracy, and even less is known about 
their evolution. 

The modified gravity approach is extremely attractive in the applications for late accelerating universe and dark 
energy. Indeed, 

1. Modified gravity provides the very natural gravitational alternative for dark energy. The cosmic speed-up is 
explained simply by the fact of the universe expansion where some sub-dominant terms (like l/R) may become 
essential at small curvature. 

2. Modified gravity presents very natural unification of the early-time inflation and late-time acceleration thanks 
to different role of gravitational terms relevant at small and at large curvature. Moreover, some models of modified 
gravity are predicted by string/M-theory considerations. 

3. It may serve as the basis for unified explanation of dark energy and dark matter. Some cosmological effects (like 
galaxies rotation curves) may be explained in frames of modified gravity. 

4. Assuming that universe is entering the phantom phase, modified gravity may naturally describe the transition 
from non-phantom phase to phantom one without necessity to introduce the exotic matter (like the scalar with wrong 
sign kinetic term or ideal fluid with EoS parameter less than —1). In addition, often the phantom phase in modified 
gravity is transient. Hence, no future Big Rip is usually expected there. 

5. Modified gravity quite naturally describes the transition from decceleration to acceleration in the universe evolution. 

6. The effective dark energy dominance may be assisted by the modification of gravity. Hence, the coincidence 
problem is solved there simply by the fact of the universe expansion. 

7. Modified gravity is expected to be useful in high energy physics (for instance, for the explanation of hierarchy 
problem or unification of GUTs with gravity). 

8. Despite quite stringent constraints from Solar System tests, there are versions of modified gravity which may be 
viable theories competing with General Relativity at current epoch. Nevertheless, more serious check of such theories 
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is necessary to fit them with various observational data and precise Solar System tests. 

These lectures are devoted to the review of some modified gravities with above mentioned attractive features and 
their cosmological properties at late universe. We do not discuss here modified gravities related with brane-world 
approach as well as Palatini formulation of modified gravities. 

The paper is organized as follows. In scct^lwe start from general f(R) gravity in accelerating FRW universe. It 
is shown how this theory may be rewritten as scalar-tensor gravity (Einstein frame). Two specific models: the one 
with l/R and i?™ terms, and another with Ini? terms are considered. It is shown that both models may lead to the 
(cosmological constant or quintessence) acceleration of the universe as well as an early time infiation. Moreover, the 
first model seems to pass the Solar System tests, i.e. it has the acceptable newtonian limit, no instabilities and no 
Brans-Dicke problem (decoupling of scalar) in scalar-tensor version. The coupling of such modified gravity with usual 
matter is considered. It is shown the occurence of transient phantom era in this case. General conditions to f{R) 
gravity admitting phantom divide crossing are derived. Finally, it is demonstrated how scalar-tensor theory may be 
rewritten as some version of modified gravity. Sect lllll is devoted to the study of FRW cosmology in string- inspired 
dilaton gravity containing scalar-Gauss-Bonnet coupling. The occurence of acceleration is again mentioned. In sect II VI 
the non-linear coupling of gravity with matter is considered. It is shown that some forms of such coupling may help 
to solve the coincidence problem simply by the fact of the universe expansion. From another point, such modified 
gravity may also predict (cosmological constant, quintessence or phantom) accelerating epoch. It is interesting that 
the model with such non-linear coupling may serve for quite interesting dynamical solution of cosmological constant 
problem. The possible applications of modified gravity for solution of hierarchy problem of high energy physics are 
described. 

In sect|3we consider /(G) and f{R,G) theories where G is Gauss-Bonnet invariant. It is shown the possibility 
of cosmic speed-up for various choices of function /, as well as possibility of decceleration-acceleration transition. 
Sect I VII is devoted to the study of late-time cosmology for the universe filled with dark fluid with inhomogeneous 
EoS. At flrst, the structure of future singularities is classified for dark fiuid with homogeneous EoS: p = —p + Ap"" . 
Next, the impact of the inhomogeneous terms to future singularities structure is analyzed for the same (effectively, 
phantom) dark fluid. As the motivation for inhomogeneous EoS dark fluid, the possible origin of the inhomogeneous 
terms from scalar-tensor or modified gravity is indicated. Some outlook and discussion are given in the last section. 



II. THE MODIFIED j(R) GRAVITY 



In the present section /(i?) gravity in FRW universe is considered. Special attention is paid to the versions of 
such theory which may describe cosmic speed-up at late universe. Let us start from the rather general 4-dimensional 
action: 

5 = y d^x^a{!(R) + L^} . (1) 

Here R is the scalar curvature, /(i?) is an arbitrary function and Lm is a matter Lagrangian density. The equation 
of the motion is given by 

- |9p./(i?) - R^.r{R) ~ V^V,/'(i?) - g^.V2/'(i?) + It^, . (2) 

With no matter and for the Ricci tensor i?^^ being covariantly constant, the equation of motion corresponding to the 
action |^ is: 

= 2/(i?) - Rf\R) , (3) 

which is the algebraic equation with respect to R. If the solution of Eq.Q is positive, it expresses deSitter universe 
and if negative, anti-deSitter universe. In the following, the metric is assumed to be in the FRW form: 

3 

ds" = -dt^+ a{tf {dx') ' . (4) 
1=1 

Here we assume that the spatial part is flat as suggested by the observation of the Cosmic Microwave Background 
(CMB) radiation. Without the matter and in FRW background, Eq.Q gives 

= -\f{R) + 3[h^ + h) f'{R) 6|/" (i?) - ISff^l ( ) /" (i?) . (5) 
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Here R is given by R = 12_ff^ + 6H. Our main purpose is to look for accelerating cosmological solutions of the 
following form: de Sitter (dS) space, where H is constant and a{t) oc e^', quintessence and phantom like cosmologies: 



when ho > (quintessence) 



ao [ts — t) ° , when ho < (phantom) 

The cosmologies of above type in various modified gravity models will be often discussed in this work. 
Introducing the auxiliary fields, A and B, one can rewrite the action |^ as follows: 



(6) 



d*3 



1 



{B{R-A) + f{A)}+C 

matter 



One is able to eliminate B, and to obtain 



d^xV^ — {f{A) (R-A)+ f{A)} + C^^tter 



(7) 



(8) 



and by using the conformal transformation g^^ ^'^ gtiu ("■ = ^ ln/'(^))) the action ((SJ is rewritten as the Einstein- 
frame action: 



1 



Se = I —2[^ - ^a'^dpad^a ~ V{a) ) + A^.t^cr 



Here, 



V{a) = e-G (c-) - e^-/ (G (e-)) = 



A 



r{A) f'iAy 



(9) 



(10) 



The action ^ is called the Jordan-frame action. In the Einstein-frame action, the matter couples with the scalar 
field a. Due to the scale transformation, there is a difference in time interval between the Einstein and Jordan frames. 
Of course, we usually use clock to measure time. The clock could be based essentially on electromagnetism. In the 
Jordan frame, where the scalar field a does not couple with matter, especially with electromagnetic fields, the time 
measured by the clock is proportional to the cosmological time in the FRW metric. On the other hand, in the Einstein 
frame where the scalar field a couples with matter, the time measured by the clock is not always proportional to the 
cosmological time, due to the coupling of the scalar field a with electromagnetic field. The difference, of course, comes 
from the factor of the scale transformation, which includes the scalar field a. Since a is time-dependent in general, 
the time measured by the clock becomes, say, slower than the cosmological time in the Einstein frame. It is estimated 
the acceleration of the universe started about five billion years ago. The five billion years are measured by using the 
electromagnetism. Hence, these five billion years could be identified with five billion years in the Jordan frame but 
not with those in the Einstein frame (for further discussion of the relation between Einstein and Jordan frames in 
cosmology, see review Q)- 

One may define the effective EoS parameter w^s in Jordan frame as 



WeS = - = -1 
P 



2H 
MP 



(11) 



The scale factor in Einstein frame (when the two frames appear) is denoted as a{t). Note that if we dont make 
transformation from one frame to another one, then scale factor will be always denoted without hat. 

A. Modified gravity with negative and positive powers of the curvature 

As the first gravitational alternative for dark energy we consider the following action Q 



f{R)^R- 



{R^AiY 



biR-A^y 



(12) 



Here we assume the coefficients n, to, c, 6 > but n, m may be fractional. Let us show that above model leads to 
acceptable cosmic speed-up and is consistent with Solar System tests. 
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For the action (|f 2|l . Eq.Q has the following form: 



= -i? + (1^^ + - 2)fe (i? - A2)™ . (13) 



Especially when n = 1 and to = 2, one gets 



If c > 0, one solution corresponds to deSitter space and another to anti-deSitter. If — ^ < c < and Ai > 0, both 
of solutions express the deSitter space. Hence, the natural possibility for the unification of early-time inflation with 
late-time acceleration appears (compare with proposal Si]). Here, higher derivative terms act in favour of early-time 
inflation while l/R term supports cosmic speed-up. 

By assuming the FRW universe metric one may define the Hubble rate by _ff = a/ a. The contribution from 
matter may be neglected. Especially when n = 1, m = 2, and Ai = A2 = in p2|) and the curvature is small, we 
obtain the following solution of 10): a cx i^, which is consistent with the result in 6]. If the present universe expands 
with above power law, the curvature of the present universe should be small compared with that of the deSitter 
universe solution in (|14l) with Ai = 0. As the Hubble rate in the present universe is (lO^'^'^eV) , the parameter c, 

which corresponds to /x^ in ji^, should be much larger than (lO~'^^eV)^. 

We now consider the more general case that f{R) is given by H12I) when the curvature is small. Neglecting the 

im ■ - (ti+i)(2,i+i) 

contribution from the matter again, solving |j2||, we obtain a (x t "+2 . It is quite remarkable that actually any 
negative power of the curvature supports the cosmic acceleration. This gives the freedom in modification of the model 
to achieve the consistency with experimental tests of newtonian gravity. 

On the other hand, when the scalar curvature R is large, one obtains a cx t ™-2 . When ^"'^^-2"^^"* 
the universe is shrinking but if we change the arrow of time by t ^ tg — t, the inflation occurs with the inverse power 
law and at t = t^, the size of the universe diverges. It is remarkable that when m is fractional (or irrational) and 
1 < TO < 2, the expression of a is still valid and the power becomes positive, the universe evolves with the (fractional) 
power law expansion. 

It is interesting that the model of this section seems to be more consistent than the one of ref. (for discussion of 
newtonian limit in theory see Q) as it may pass Solar System tests. Indeed, in small gravitational object like 
the Earth or the Sun in the model [3 is considered. It has been shown that the system becomes instable. This may 
cause the unacceptable force between distant Galaxies 0|. As shown in 0, however, by adding the positive power 
(higher than 1) of the scalar curvature to the action, the instability could be significally improved. Moreover, the 
account of quantum effects in modified gravity also acts against instability 

It has been also mentioned in that l/R model jj^l which is equivalent to some scalar-tensor gravity is ruled 
out as realistic theory due to the constraints to such Brans-Dicke type theories. In iJl, it was shown that by adding 
the scalar curvature squared term to the action, the mass of the scalar field can be adjusted to be very large and 
the scalar field can decouple. Hence, the modified gravity theory of this section (after some fine-tuning passes 
the Solar System tests. In addition, at precisely the same values of parameters the above modified gravity has the 
newtonian limit which does not deviate significally from the one in General Relativity. 



B. InR gravity 

Other gravitational alternatives for dark energy may be suggested along the same line. As an extension of the 
theory of the previous section, one may consider the model 12] containing the logarithm of the scalar curvature R: 

f{R) = R + a'\n^+f3R"' . (15) 
M 

We should note that m = 2 choice simplifies the model. Assuming the scalar curvature is constant and the Ricci 
tensor is also covariantly constant, the equations ||2Jl are: 

R 

= 2/(i?) - RJ'{R) = j{R) EE i? + 2a' In — - a' . (16) 

M 

If a' > 0, f{R) is monotonically increasing function and lim^^o f{R) = —00 and lim/f^_|_oo f{R) = +00. There is 
one and only one solution of (|16|) . This solution may correspond to the inflation. On the other hand, if a' < 0, 
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limfl^o f{R) = oo and liinj^_+oo = +00. Since /'(i?) = 1 + the minimum of f{R), where f'{R) = 0, is given by 
R = -2a'. If f{-2a') > 0, there is no solution of If /(-2a') = 0, there is only one solution and if /(-2a') < 0, 
there are two solutions. Since /(—2a') — —2a' (^l — In j , there are two solutions if — ^ > e. Since the square 

root of the curvature R corresponds to the rate of the expansion of the universe, the larger solution in two solutions 
might correspond to the inflation in the early universe and the smaller one to the present accelerating universe. 

We can consider late FRW cosmology when the scalar curvature R is small. Solving (0), it follows that the power 
law inflation could occur: a oc . Since a > but a < 0, the deccelerated expansion occurs. 

One may discuss further generalizations like jlJl (see also [T^ 'l 

/(i?) = i? + ^i?-«(^ln-j . (17) 

Here n is rectricted by n > — 1 (m is an arbitrary) in order that the second term could be more dominant than the 
Einstein term when R is small. For this model, we find 

(n + l)(2,i+l) 

a ~ t . (18) 

This does not depend on m. The logarithmic factor is almost irrelevant. The effective WoS is given by 

6n^ + 7n — 1 

^^'^ = -3(n + l)(2n+l) ' ^^^^ 



Then Wcs can be negative if 



1 -7+V73 , , 

-l<n<--orn> ^ — = 0.1287---. (20) 



From (fTH|l . the condition that the universe could accelerate is > ^j^^t is: 



n+2 

n> '^^^ = 0.366 ■■■ . (21) 
Clearly, the effective EoS parameter w may be within the existing bounds. 

C. Modified gravity coupled with matter 

In this subsection we will show that the modified gravity may lead to the effective phantom or quintessence dark 
energy without necessity to introduce the (negative kinetic energy) phantom scalar or negative pressure ideal fluid. 
In fact, the matter is taken to be the usual ideal fluid. 

The ideal fluid is taken as the matter with the constant w: p — wp. Then from the energy conservation law it 
follows p — pofl ■ In a some limit, strong cuvature or weak one, f{R) rnay behave as f{R) ^ foR"", with 

constant /o and a. An exact solution of the equation of motion is found to be |14| 

2a 



3(1 + w) 



PO 



6ho + 12/1^) {(1 - 2a) (1 - a) - (2 - a)ho} 



(22) 



When a = 1, the result /iq = 3(Urw) Einstein gravity is reproduced. The effective WeS may be defined by 

ho — 3(1^^77-^ • By using (|^ . one finds the effective WeS lfTT|l is given by 

Woff = -1 + ^-t^ . (23) 
a 

Hence, if w is greater than —1 (effective quintessence or even usual ideal fluid with positive w), when a is negative, we 
obtain the effective phantom phase where is less than — 1 . This is different from the case of pure modified gravity. 
Moreover, when a > w + 1 (it can be even positive), WeS could be negative (for negative w). Hence, it follows that 
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modified gravity minimally coupled with usual (or quintessence) matter may reproduce quintessence (or phantom) 
evolution phase for dark energy universe in an easier way than without such coupling. 
One may now take /(J?) as 

/(i?) = 4(^-^^'" + ^^') • (24) 

AC 

When the cuvature is small, the second term becomes dominant and one may identify fo — and a = —n. Then 
from H23|l . it follows WcS = — 1 — i±2£. Hence, if n > 0, an effective phantom era occurs even if w > —1. Usually 
the phantom era generates the future Big Rip singularity J^] . (Note that there may occur different types of future 
sing ularities, say, sudden ones etc. The complete classification of Big Rip singularities is given in 

|l8l|.'l However, near the Big Rip singularity, the curvature becomes large and the last term becomes dominant. In 
this case a = 2 and Wcs = — 'j-^. Then if w > —1, it follows Wch > —1, which prevents the Big Rip singularity 
(it makes the phantom phase transient). To conclude, it looks quite promising that modest modification of General 
Relativity coupled to ideal fluid matter leads to an effective dark energy universe in the very natural way. 

D. The equivalence with scalar-tensor theory 

It is very interesting that f{R) gravity is in some sense equivalent to the scalar-tensor theory with the action: 

6.(0) = -L/i'(0) , y(0) ^ L (3/,(0)2 + /,'(0)) . (25) 

Here /i((/)) is a proper function of the scalar field (p. Imagine the following FRW cosmology is constructed: 

= i , H = h{t) . (26) 

Then , as it has been demonstrated in ref.|23|, any cosmology defined by iJ = h{t) in (|26|) can be realized by H25|l . 
Indeed, if one defines a new field (p as 



^ = / d^VR^ , (27) 



the action H25|l can be rewritten as 

S = J d'x^ { T \d^vd^v -V{^)^ . (28) 

The sign in front of the kinetic term depends on the sign of w(0). If the sign of u and therefore the sign of H is 
positive (negative), the sign of the kinetic term is — (+). Therefore, in the phantom phase, the sign is always + and 
in the non-phantom phase, always — . One assumes can be solved with respect to (/?: 4> = 4>{'p). Then the potential 
V(lp) is given by V{^p) = V {4>{ip)). Since V{^) could be uniquely determined, there is one-to-one correspondence 
between H and V{tf). 

In case the sign in front of the kinetic term of Lp in H28|l is — , we can use the conformal transformation g^^ — > 
e^'^'^v'^g^j/, and make the kinetic term of ip vanish. Hence, one obtains 



S= l__Ii?-e±^«-V§f(^) . (29) 



The action H29|) may be called as Jordan frame action and the action (|28|) as the Einstein frame action. Since ip 
becomes the auxiliary field, one may delete p) by using an equation of motion: 



(30) 



Modified Gravity 



7 



which may be solved with respect to R as ip = (p{R). One can rewrite the action 129(1 in the form of f{R) gravity |20| : 

S ^ J d^xV^fiR) , 

f{R) ^ ^——^ i?-e±2-'p(fl)V§F(^(i?)) . (31) 

2k^ 

Note that one can rewrite the scalar-tensor theory H25|l or equivalently ((28|l , only when the sign in front of the kinetic 
term is — in (|28|l . that is, is positive. In the phantom phase, uj{(l)) is negative. In this case, the above trick to 
rewrite the phantom scalar-tensor theory as f{R) gravity does not work. We should note, however, as the metric is 
conformally transformed, even if the universe in the original metric (corresponding to the Einstein frame) is in the 
phantom phase, the universe in the rescaled metric (corresponding to the Jordan frame) can be in non-phantom phase 
in general. In addition, one should bear in mind that the same FRW dynamics (same scale factor) may be obtained 
from one of these two mathematically equivalent descriptions. Nevertheless, the properties of such FRW cosmology in 
both models may be quite different (for instance, more FRW solutions may exist, say, for scalar-tensor theory, newton 
law may be different, stability of such equivalent solutions may be different, etc). 



E. Transition from phantom to non-phantom phase 



Assuming that our universe is entering into the phantom phase recently, and having in mind the occurence of (cos- 
mologically recent) decceleration-acceleration transition, the question appears: what are the conditions for phantom- 
non-phantom transition in modified /(i?) gravity? For instance, the conditions for the existence of deSitter universe 
solution in J{R) gravity may be derived in the form of algebraic equation '21^. The stability of deSitter universe 
solution was established in ref.|2l| where one-loop quantization of f{R) was done. (The same stability criteria was 
also found in ref.j^^l using perturbations analysis.) 

The effective EoS parameter is defined by (|ll|l . We now define a time-dependent parameter F{t) by 

no^-g,. (32) 

that is, F{t) = 1-1- Wcff. Note that = at the transition point between non-phantom (wcff > — 1 or F > 0) and 
phantom (wes < — 1 or _F < 0) phases. 

We now consider the region near the transition point between non-phantom and phantom phases and assume F is 
small. Then Eq.Q reduces to 

= -i/(12i/2) + 3H^f{12H^) + 81H'^Ff"{12H^) + 27H^Ff" {12H^) + O (F^) . (33) 

If the transition point corresponds to the deSitter solution, which satisfies 

= -i/(12ij2) + 3H^f'{12H^) , (34) 

one finds F ^ ~3HF. Then F behaves as ~ FqC"'^^*, which does not vanish in a finite time and therefore there 
does not occur the transition. On the other hand, if the transition point does NOT correspond to the deSitter solution, 
it follows 

• f{12H^)-6H^f{12H^) 

54i73/"(12i72) • ^'^''> 



Then near the transition point if 



/"(12i?2) > and f{l2H^) - QH^f'{\2H^) < 



or f"{l2H') < and /(12i/") - QH^f'{l2H^) > , (36) 

there could occur the transition from the non-phantom phase {F > 0) to phantom phase (F < 0). 
As an example, we consider 

fiR) = /oe^/«^o , (37) 
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with constants /q and Hq (for similar choice with negative R in the exponent see H^.) From (|35|l . one finds 

F '-■ (38) 

Then by properly choosing initial conditions, \i H > Hq when F ^ but positive (that is, in non-phantom phase), 
we find F < and there could occur the transition from non-phantom phase to phantom phase (crossing of phantom 
divide). 

III. STRING-INSPIRED GAUSS-BONNET GRAVITY AS DARK ENERGY 

The successful dark energy theory may be searched in string/M-theory. Indeed, it is quite possible that some un- 
usual gravity-matter couplings predicted by the fundamental theory may become important at current, low-curvature 
universe (being not essential in intermediate epoch from strong to low curvature). For instance, in the study of string- 
induced gravity near to initial singularity the role of Gauss-Bonnet (GB) coupling with scalar was quite important 
for ocurrence of non-singular cosmology |2J,|23. Moreover, string/M-theory naturally predicts the appearence of the 
terms with inverse powers of curvature invariants in compactified, low-energy effective action l2^. 

The present section is devoted to the study of the role of GB coupling with the scalar field (which naturally appears 
in low-energy effective action) to the late-time universe ■ It is explicitly demonstrated that such term itself cannot 
induce the effective phantom late-time universe if the scalar is canonical in the absence of potential term. It may 
produce the effective quintessence (or phantom) era, explaining the current acceleration only when the scalar is 
phantom or when the scalar is canonical with non-zero potential. It is interesting that it may also have the important 
impact to the Big Rip singularity, similarly to quantum effects 0,0,123, preventing it. Note that we concentrate 
mainly on the exponential scalar-GB coupling and exponential scalar potential, while the consideration of other types 
of such functions and their role in late time cosmology may be considered too p9| . 

We consider a model of the scalar field coupled with gravity. As a stringy correction, the term proportional to 
the GB invariant G = — AR/^^R^^" + R^.^paR^"'"^ is added. The starting action is given by 

V = Voe-^ , /(0) = foe^ . (39) 

Here 7 — ±1. For the canonical scalar, 7=1 but at least when GB term is not included, the scalar behaves as 
phantom only when 7 = — 1. 

Starting with FRW universe metric Ullfl in scct^ and assuming JHl in sectEI the following solutions may be 
obtained 

^o.?^-.^(3.g(l-M + ^^^^^^}, 

(40) 



(1 + ho) 

48fohl _ 6 / j(f>lK 



q k2 (1 + ho) 

Even if 7 = —1, there appear the solutions describing non- phantom cosmology coresponding the quintessence or 
matter. As an example, we consider the case that Hq = — ^ < —1, which gives Wcff = —1.025. Simple tuning gives 
other acceptable values of the effective w in the range close to —1. This is consistent with the observational bounds 
for effective w (for recent discussion and complete list of refs., see H^])- Then from H40|l . one obtains 

2 1 /531200 403 , 2\ fo I f 9 27 , 2\ 

Vot^ = h 70oK^ , ^ = ^ H -fd}on ■ 41 

" ^ k2 1^ 231 154 '™ J tj k2 V49280 7884800 '™ ^ ^ ' 

Therefore even starting from the canonical scalar theory with positive potential, we may obtain a solution which 
reproduces the observed value of w. 

There is another kind of solutions j23|. If ^ and H are constants: (j) ^ tpo, H = Hq, this corresponds to deSitter 
space. Then the solution of equations of motion gives: 
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Therefore in order for the solution to exist, the condition is /o < 0. In 1)42(1 . (po can be arbitrary. Hence, the Hubble 
rate H — Hq might be determined by an initial condition. 

In 27] , the stability of the above (and related) solutions was checked and it was found that the case corresponding 
to phantom cosmology with /iq < is always stable but the case corresponding to non-phantom cosmology with Hq > 
is always unstable. Also in 31J, by taking into account the higher-order string corrections to Einstein-Hilbert action, 
the evolution of (phantom) dark energy universe was studied. When fixed dilaton, while the presence of a cosmological 
constant gives stable de-Sitter fixed points in the cases of heterotic and bosonic strings, no stable de-Sitter solutions 
exist when a phantom fluid is present. It was found that the universe can exhibit a Big Crunch singularity with a finite 
time for type II string, whereas it reaches a Big Rip singularity for heterotic and bosonic strings. Thus the fate of 
dark energy universe crucially depends upon the type of string theory under consideration. Furthermore, a barotropic 
perfect fluid coupled to the scalar field (dilaton or compactification modulus) was also considered and phase space 
analysis and the stability of asymptotic solutions are performed for a number of models which include (j) fixed scalar 
field, (a) linear dilaton in string frame, and (iii) logarithmic modulus in Einstein frame. It was also shown that 
dilaton Gauss-Bonnet gravity (with no matter) may not explain the current acceleration of the universe. It was also 
studied the future evolution of the universe using the GB parametrization and found that Big Rip singularities can be 
avoided even in the presence of a phantom fiuid because of the balance between the fiuid and curvature corrections. 
A non-minimal coupling between the fluid and the modulus fleld also opens up the interesting possibility to avoid Big 
Rip regardless of the details of the fluid equation of state. Note also that above action has interesting applications in 
early time string cosmology (for recent discussion and list of references, see 32] ). Definitely, dilaton Gauss-Bonnet 
gravity (perhaps, also with account of higher order terms like Euler invariant) deserves very deep investigation as 
possible candidate for dark energy. 

IV. MODIFIED GRAVITY: NON-LINEAR COUPLING, COSMIC ACCELERATION AND 

HIERARCHY PROBLEM 



It is interesting to investigate the role of non-minimal coupling of modified gravity with dark ener gy to cosmic 
acceleration epoch. The corresponding, quite general theory has been suggested in ref.|33j (see also |3J])- an 
example of such theory, the following action is considered: 



Here Ld is matter-like action (dark energy). The choice of parameter /i may keep away the unwanted instabilities 
which often occur in higher derivative theories. The second term in above action describes the non- linear coupling 
of matter with gravity (in parallel with R(/)^ term which is usually required by renormalizability condition 35]). 
Similarly, such term may be induced by quantum effects as non-local effective action. Hence, it is natural to consider 
that it belongs to matter sector. (Standard matter is not included for simplicity). 
By the variation over y^j,, the equation of motion follows: 



A. Gravitational solution of coincidence problem 




(43) 




(44) 



Here the effective EMT tensor T^^ is defined by 





(45) 



Let free massless scalar be a matter 



(46) 



Then the equation given by the variation over (/) has the following form: 




(47) 



Modified Gravity 



10 



The metric again corresponds to FRW universe with flat 3-space. If we assume </) only depends on i (0 = 4'{t)), the 
solution of scalar field equation 147|l is given by 

= qa-^R-° . (48) 

2 

Here q is a constant of the integration. Hence R"Ld — 2^^, which becomes dominant when R is small (large) 
compared with the Einstein term -^R ifo; > —1 (a < —1). Thus, one arrives at the remarkable possibility 'SS] that 
dark energy grows to asymptotic dominance over the usual matter with decrease of the curvature. At current universe, 
this solves the coincidence problem (the equality of the energy density for dark energy and for matter) simply by the 
fact of the universe expansion. 

Substituting H48(l into ()44|l . the (/x, ^) = (i,t) component of equation of motion has the following form: 



ft. 



(^6H+12H^ 



13 2\ A„2 7 



—a + a'j HH' + (^1 + ^ J • (49) 

The accelerated FRW solution of exists 

f„ " + 6 _ K^q' (2a - 1) (g - 1) 
a = ant 3 [ = — 7, — , On = -; . (50) 

V 3t ; ' ° ;,2o3(«+l)"+l(2(2^_l))-+2 ^ ' 

Ea. (|5()|l tells that the universe accelerates, that is, a > if a > 2. If a < —1, the solution H50|l describes shrinking 
universe if t > 0. If the time is shifted as t t ~ {t^ is a constant), the accelerating and expanding universe 
occurs when t < tg. In the solution with a < — 1 there appears a Big Rip singularity at i = tg. For the matter 
with the relation p = wp, where p is the pressure and p is the energy density, from the usual FRW equation, one has 
a oc i . For a oc t'^" it follows w = —1 + g^^, and the accelerating expansion {ho > 1) of the universe occurs if 
— l<w< — i. For the case of (O, one finds 

(51) 



1 



Then if a < —1, we have w < —I, which is an effective phantom. For the general matter with the relation p ~ wp 
with constant w, the energy E and the energy density p behave as E ^ a~^^ and p ~ a~^^'^~^^\ Thus, for the 
standard phantom with w < —1, the density becomes large with time and might generate the Big Rip. Hence, 
it is demonstrated that non-linear coupling of gravity with matter may produce the accelerated cosmology with the 
effective EoS parameter being close (above, equal or below) to —1. More complicated couplings 's^ may be considered 
in the same fashion. 



B. Dynamical cosmological constant theory: an exact example 



There are many proposals to solve the cosmological constant problem dynamically. One of the possible solutions 
of the cosmological constant problem is pointed out by Mukohyama and Randall ,3^ (see also [^). In [s^, the 
following action similar to the one under consideration has been proposed: 

where /(i?) is a proper function. When the curvature is small, it is assumed f{R) behaves as 

/(i?) ^ {.^rT ■ (53) 

Here m is positive. When the curvature is small, the vacuum energy, and therefore the value of the potential becomes 
small. Then one may assume, for the small curvature, V{ip) behaves as 

Vi^p) ^Voip- Vc) . (54) 



R 

2^ 



aoR^ 



2qK^f{RY'i- 



(52) 
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Here Vq and ipc are constants, li q > 1/2, the factor in front of the kinetic term of (p in (|52|l becomes large. This 
makes the time development of the scalar field ip very slow and it is expected that ip does not reach ipc- This model 
may explain the acceleration of the present universe. The model H52|l is also expected to be stable under the radiative 
corrections. In fact, when the curvature is small and the time development of the curvature can be neglected, if we 
rescale the scalar field curvature in the kinetic term can be absorbed into the redefinition of 

p and there appear factors including Therefore if m{2q — l)/q > 0, the interactions could be suppressed 

when the curvature is small and there will not appear the radiative correction to the vacuum energy except the 
one-loop corrections. 

There is an exactly solvable model which realizes the above scenario. Let us choose 

f{R) = /3i?2 , V{p) = Vo{p- Pc) . (55) 

Here /3 is a constant. term is neglected by putting = in H52|) since the curvature is small. Searching for the 
solution ^ in sect Inland choosing p ^ pc + po/t'^ or p — pc + po / {ts ~ t)'^ , the following restrictions are obtained 

m 

5Ap{-l + 2hofht 3/10 + 1 

^0 «2 ^i2hl - 2h, - 1) ' ^6K^12hl-2ho-l){-l + 2ho) ' ^ ^ 

Since Pq should be positive, one finds 

u ^ n l-\/l3 , 1 + VT3 ,^1 

when /3 > , — — < ho < — — or ho > - , 

when /3 < , ho < ^ or ^ < ho < ^ . (57) 

For example, if ho = —1/60, which gives Wcs = —1.025, we find 

K^^o = ±^y^ = ±0.388722... . (58) 

It is interesting that the value of WcS is consistent with the observed bounds. (Note that recently the above model 
has been reconsidered in ref. 39] ). 

For ho > case, since R = 6H + 12H^, the curvature R decreases as t^^ with time t and p approaches to pc but 
does not arrive at pc in a finite time, as expected in 's^. 

As H behaves as ho/t or ho/{ts — t) for © in scctEl if we substitute the value of the age of the present universe 
10^°years~ (10~^'^eV)~^ into t or tg — t, the observed value of H could be reproduced, which could explain the 
smallness of the effective cosmological constant A ~ H^. Note that even if there is no potential term, that is, Vq = 0, 
when /3 < 0, there is a solution 

/lo - -J <^^ = -0.2171... , (59) 

which gives the EoS parameter : w — —3, although w is not realistic. Playing with different choices of the potential 
and non-linear coupling more realistic predictions may be obtained. 



C. Hierarchy problem in modified gravity 

Hierarchy problem is known to be the fundamental one in high energy physics. It is interesting to understand if 
modified gravity approach may help in resolution of hierarchy problem? Recently the hierarchy problem has been 
investigated, in j40], in frames of scalar-tensor theory. Here, similar proposal will be discussed in the modified gravity. 
In (40j |. in order to generate the hierarchy, a small scale, which is the vacuum decay rate Fuac, was considered. Instead 
of r„ac5 we here use the age of the universe ~ 10"^'^ eV, as the small mass scale. 

We now start from the action for modified gravity coupled with matter (Q) in sectm Let us assume that the matter 
Lagrangian density £mattcr contains a Higgs-like scalar field p 

1 V? 

/:matter = -2V^<^V^<^+ Y'/''- V + -- - ■ (60) 
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Under the conformal transformation g^^ — !■ e'^ g^u, the matter Lagrangian density £mattcr is transformed as 



^matter ^ rSiattcr = -^V^(^V^(^ + ^^^^ _ ^^2.^4 ^ . . . _ (g^) 



2 ^ ■ 2 

By redefining as — > oT" /"^Lp^ C^^^^^^ acquires the following form 



1 V?(f 



^matter ~ "TT^VV^^ + - A^* + • • • , (62) 



where the time derivative of cr is neglected. Then, the massive parameter /i, which determines the weak scale, is 
effectively transformed as /i ^ /i = e'^/^/i. In principle, /i can be of the order of the Planck scale 10"'^^ GeV, but if 
~ 10~^^ in the present universe, \1 could be of the order of 10^ GeV, which is the scale of the weak interaction. 
Therefore there is a quite natural possibility that the hierarchy problem can be solved by above version of modified 
gravity (similar solution of hierarchy problem may be found in two scalar-tensor gravity of ref.|4l| as it has been 
shown in 42].) 

One may consider the model (42| : 

/(i?) = i? + /oi?" , (63) 

with constant /q and a. (Without the first term the restrictions to the parameter a have been recently studied in 
ref.|4.l|.) If a < 1, the second term dominates, when the curvature is small. Assuming that the EoS parameter, w, of 
matter is constant one gets 14] the accelerated FRW cosmology solution with the effective EoS parameter Wcff (|^ . 
Neglecting the first term in H63() . it follows 

e"/2 ^ 1 . (64) 

In the present universe, R - (10^33 eV) . Assume now that /q is of the order of the Planck scale - 10^^ GeV = 10^** eV 
as /o ~ (lO^^eV)^*^. Then, Eq. ^ yields 

^ io6i("-i) . (65) 

If we furthermore assume that —17 = 61 (a — 1), we find a = 44/61. In that case, if w > —1, Wcff > —1 and 
the universe is not in its phantom phase, but (|65|l hints towards the possibility that modified gravity can solve the 
hierarchy problem. 

Let us write the action of the scalar-tensor theory as 

S ^ c?'a;^/^e"'^(i?- ia^</.9'^</)-yoe-'*/^°)) + J d^x^C^,,,,, , 

,2 



i5Va^^-t-^(p2-A^^ + --. . (66) 



This action could be regarded as the Jordan frame action (see sectnj- Scalar field could be identified with the 
Higgs field in the weak(-electromagnetic) interaction. Now the ratio of the inverse of the the effective gravitational 
coupling k — Ke~~^ and the Higgs mass fj, is given by 

1 

^ C 2 

^ = . (67) 

Hence, even if both of 1/k and ^ are of the order of the weak interaction scale, if e^ 10"'^^, 1/k could have an order 
of the Planck scale. 

By rescaling the metric and the Higgs scalar (p as 

g^.^ ^ e-°"*'gf,^ , ip^e'^ip, (68) 

the Einstein frame action is obtained: 



d X-\J .9-^ matter 5 

+ \d^<i>v) (d^^ + \dAv) + - + • • • . (69) 
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In the Einstein frame, the gravitational couphng k is constant but the effective Higgs mass jl defined by /2 = e 2 /i 
can be time-dependent. Hence, the ratio of 1/k and /2 is given by 

1 

4 - — , (70) 

which is identical with (|67|l . Hence, even if both of 1/k and fx are of the order of the Planck scale, if e^ ^ 10^^, jl 
could have an order of the weak interaction scale. Therefore the solution of the hierarchy problem does not essentially 
depend on the frame choice. 

Nevertheless, note that the cosmological time variables in two frames could be different due to the scale transfor- 
mation (|68|l as dt ^ di = c^~dt. Therefore the time- intervals are different in the two frames. The units of time and 
the length are now defined by the electromagnetism. Then the frame where the electromagnetic fields do not couple 
with the scalar field ((> could be physically more preferrable. Since the electromagnetic interaction is a part of the 
electro-weak interaction, the Jordan frame in (|66|l could be more preferrable from the point of view of the solution of 
hierarchy problem. 

V. LATE-TIME COSMOLOGY IN MODIFIED GAUSS-BONNET GRAVITY 

A. /(G) gravity 



Our next example is modified Gauss-Bonnet gravity. Let us start from the action |42, |4J| : 



S - / d^xV^ ( + f{G) + Cm ] . (71) 



2^2 

Here Cm is the matter Lagrangian density and G is the GB invariant: G = ~ AR^^^R^" + R^^^^^R'^'^^'^ . By variation 
over g^^ one gets: 

= -i^ (-R"'' + lg'"'R] + T'"' + Ig'^'fiG) - 2f'(G)RR'"' 
\ 2' J 2' 

+Af'{G)R''pR''P - 2f'{G)R^''"'^R%„^ - 4f'{G)R''P''''Rp„ + 2 {V^V^fiG)) R 
-2(7'"' (V2/'(G)) i? - 4 {VpVfiG)) R^P - 4 (VpV''/'(G)) ff"' 

+4 (V V(G)) i?^"' + 45^'' (VpV^/'(G)) RP" - 4 (VpV„/'(G)) R'^p'"' , (72) 

where T^''^ is the matter EM tensor. By choosing the spatially-flat FRW universe metric (@J) in scctm the equation 
corresponding to the first FRW equation has the following form: 

= -\h^ + Gf'{G) - ,f{G) - 2AGf"{G)H^ + pm , (73) 

where pm is the matter energy density. When pm — 0, Eq. 173|l has a deSitter universe solution where H , and therefore 
G, are constant. For H = Hq, with constant Hq, Eq. (|73|l turns into 

= - + 24i7„V' (24i/o') - / (24i?o') ■ (74) 

For a large number of choices of the function f{G), Eq. (|74ll has a non-trivial {Hq ^ 0) real solution for Hq (deSitter 
universe). The late-time cosmology for above theory without matter has been discussed for a number of examples in 
refs. [4^ . We follow these works here. 

We now consider the case pm 7^ 0. Assuming that the EoS parameter w = Pm/ Pm for matter {pm is the pressure 
of matter) is a constant then, by using the conservation of energy: pm + 3iJ {pm + Pm) — 0, we find p = poa~^'^^~^^\ 
The function /(G) is chosen as 

/(G) = /o|G|'' , (75) 

with constant /o and p. If /3 < 1/2, /(G) term becomes dominant compared with the Einstein term when the 
curvature is small. If we neglect the contribution from the Einstein term in (j73|l . assuming FRW anzats ((HJ in sectHll 
the following solution may be found 



"•0 — TTTi — ; — \ ' ^0 



3(1 -fH ' " L {ho~l)po 



/o(/3-l) . , . Ml/3 



{2'i\hli-l + ho)\Y {ho~l + A(3) . (76) 
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Then the effective EoS parameter Wcs (|ll|l in sect ^ is less than — 1 if /? < 0, and for it; > —1 is 

which is again less than —1 for /3 < 0. Thus, if < 0, we obtain an effective phantom with negative ho even in the 
case when w > —1. In the phantom phase, there might seem to occur the Big Rip at i = 15] . Near this Big Rip 
, however, the curvature becomes dominant and then the Einstein term dominates, so that the /(G) term can be 
neglected. Therefore, the universe behaves as a = oot^/'^*-'""'"^-' and as a consequence the Big Rip does not eventually 
occur. The phantom era is transient. 

The case when < /? < 1/2 may be also considered. As (3 is positive, the universe does not reach here the phantom 
phase. When the curvature is strong, the /(G) term in the action H71|l can be neglected and we can work with 
Einstein gravity. Then if w is positive, the matter energy density should behave as Pm ^ t^'^, but /(G) goes 
as /(G) ~ t~'^^. Then, for late times (large t), the /(G) term may become dominant as compared with the matter 
one. If we neglect the contribution from matter, Eq. H73|) has a deSitter universe solution where H, and therefore 
G, are constant. If _ff = Hq with constant Hq, Eq. H73() looks as H74() . As a consequence, even if we start from the 
deceleration phase with w > —1/3, we may also reach an asymptotically deSitter universe, which is an accelerated 
universe. Correspondingly, also here there could be a transition from acceleration to deceleration of the universe. 

Now, we consider the case when the contributions coming from the Einstein and matter terms can be neglected. 
Then, Eq. reduces to 

= G/'(G) - /(G) - 24Gf'{G)H^ . (78) 
If /(G) behaves as (I75II . from assumption JB)) in scct^ we obtain 

- (/? - 1) hi {ho - 1) (/lo - 1 + 4/?) . (79) 



As ft-o = 1 implies G — 0, one may choose 



and Eq. lfTT)l in sect HTl gives 



/lo = 1 - 4/3 , (80) 



-e. = -l + ^^^. (81) 

Therefore, if /3 > 0, the universe is accelerating {wos < —1/3) and if /3 > 1/4, the universe is in a phantom phase 
{wcs < —1). Thus, we are led to consider the following model: 

/(G) = /. \Gf' + fi \Gf' , (82) 

where we assume that 

A > ^ , ^ > > i ■ (83) 

Then, when the curvature is large, as in the primordial universe, the first term dominates, compared with the second 
one and the Einstein term, and gives 

-l>z.e. = -l + ^(^>-5/3. (84) 

On the other hand, when the curvature is small, as is the case in the present universe, the second term in l|82f) 
dominates, compared with the first one and the Einstein term, and yields 

^e.^-l + ^^^<-5/3. (85) 

Therefore, theory 1)82(1 can in fact produce a model which is able to describe inflation and the late-time acceleration 
of the universe in the unified way. 

Instead of H83|) . one may also choose /3; as 

l>Pi>0, (86) 
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which gives 

> Woff > -1 . (87) 

Then, what we obtain is an effective quintessence. Moreover, by properly adjusting the couplings fi and // in H82|l . 
one can obtain a period where the Einstein term dominates and the universe is in a deceleration phase. After that, 
there would come a transition from deceleration to acceleration, when the GB term becomes the dominant one. More 
choices of f{G) may be studied with the purpose of the construction of natural accelerating universes. 

Let us address the issue of the correction to Newton law. Let g(o) be a solution of H72() and represent the perturbation 
of the metric as g^^, = g[Q)^^ + h^^^. First, we consider the perturbation around the deSitter background which is a 
solution of H74(l . We write the deSitter space metric as g(o)fiv, which gives the following Riemann tensor: 

R(id)tiiypa = Hi {g(Q)ppg(Q)ycr - 9 (0) tia 9 {Q) v p) ■ (88) 

The fiat background corresponds to the limit of Hq — > 0. For simplicity, the following gauge condition is chosen: 
9(o)'^t^u = V^o)/i^^ = 0. Then Eq. ^ gives 



= ^ (V^/i^, - 2Hlh^,) + T^, . (89) 



The GB term contribution does not appear except in the length parameter I/Hq of the deSitter space, which is 
determined with account to the GB term. This may occur due to the special structure of GB invariant. Eq. H89|l tells 
us that there is no correction to Newton's law in deSitter and even in the flat background corresponding to Hq 0, 
whatever is the form of / (at least, with above gauge condition). Note that the study of newtonian limit in 1/R 
gravity (where significant corrections to newton law may appear) and its extensions has been done in 0, Q . For most 
1/R models the corrections to Newton law do not comply with solar system tests. 
By introducing two auxilliary fields, A and B, one can rewrite the action (|71|l as 



S = 



J d^^V^{^R + B{G^A) 

+/(A) + /:„) . (90) 

Varying over B, it follows that A = G. Using this in (|90|l . the action H71|l is recovered. On the other hand, varying 
over A in one gets B = f'{A), and hence 

d'xV^[^R + f'{A)G - AfiA) + fiA)) . (91) 



By varying over A, the relation A = G is obtained again. The scalar is not dynamical as it has no kinetic term. We 
may add, however, a kinetic term to the action by hand 



S- 



d^Ad^A + f'iA)G - AfiA) + /(A)) . (92) 



Then one obtains a dynamical scalar theory coupled with the Gauss-Bonnet invariant and with a potential. It is 
known that a theory of this kind has no ghosts and it is stable, in general. Actually, it is related with string-inspired 
dilaton gravity proposed as alternative for dark energy |23,|21|- Then, in case that the limit e — > can be obtained 
smoothly, the corresponding f{G) theory has no ghost and could actually be stable. It may be of interest to study 
the cosmology of such theory in the limit e ^ 0. 



B. /(G, i?) gravity 

It is interesting to study late-time cosmology in generalized theories, which include both the functional dependence 
from curvature as well as from the Gauss-Bonnet term |43 |: 

S= f d^x^jifiG,R) + C,n) ■ (93) 
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The following solvable model is considered: 



The FRW solution may be found again: 

ho 



^ - 2/0 ± Jsfo /o - ^) 
^ , ho = ^ — . (95) 



Then, for example, if K^fo < —3, there is a solution describing a phantom with < —1 — \/2 and a solution describing 
the effective matter with ho > — 1 + \/2. Late-time cosmology in other versions of such theory may be constructed. 



VI. INHOMOGENEOUS EQUATION OF STATE OF THE UNIVERSE DARK FLUID 

Let us remind several simple facts about the universe filled with ideal fluid. By using the energy conservation law 
— p + 3H {p + p), when p and p satisfy the following simple EOS p — wp with constant w, we find p = poa~^'^^^^\ 

Then by using the first FRW equation {3/k^)H^ — p, the well-known solution follows a — ao {t — ti) ^("^'+1' {w > — 1) 

or ao {t2 — t) w ^ —1 (w < —1) and a — aoe^'V^ when w = —1, which is the deSitter universe. Here ti and t2 

are constants of the integration. When w < —I, there appears a Big Rip singularity in a finite time at t = ^2- 

In general, the singularities in dark energy universe may behave in a different way. One may classify the future 
singularities as follows 0: 

• Type I ("Big Rip") : For t ^ ts, a 00, p ^ 00 and \p\ 00 

• Type II ("sudden") : For t ^ ts, a as, p — > Ps or and \p\ — )■ 00 

• Type III : For t ~f tg, a— >as, p 00 and \p\ — > 00 

• Type IV : For t ^ tg, a ^ as, p ^ 0, |p| ^ and higher derivatives of H diverge. This also includes the case 
when p (p) or both of them tend to some finite values while higher derivatives of H diverge. 

Here ts, fls and ps are constants with 7^ 0. The type I may correspond to the Big Rip singularity, which emerges 
when constant w is less than —1. The type II corresponds to the sudden future singularity 16] at which a and p are 
finite but p diverges. The type HI appears for the model with p = —p — Ap'^ J^, ,2SJ, which is different from the 
sudden future singularity in the sense that p diverges. Now, let us discuss the universe filled with ideal fluid with 
inhomogeneous EoS. (It is assumed that GR is valid). Effectively, this approach overlaps with modified gravity one as 
inhomogeneous contribution may be considered as one coming from new gravitational terms. Moreover, it may be an 
extremely useful in the situation when modified gravity represents highly non-linear theory where only some features 
(but not the whole Lagrangian) are known like is currently the case with M-theory. 

A. The singularities in the inhomogeneous EoS dark fluid universe 

One may start from the dark fluid with the following EOS: 

P--P-/(P), (96) 

where /(p) can be an arbitrary function in general. The choice /(p) oc p" with a constant a was proposed in R.ef. |28|| 
and was investigated in detail in Ref . J^] . Then the scale factor is given by 



1 ( dp_\ 



and the cosmological time may be found 



"""'""[i^J m) • 



As an example we may consider the case that 

/(p) = Ap" . (99) 

Then we find El: 
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• In case Q; = l/2ora = 0, there does not appear any singularity. 

• In case a > 1, when t —i- Iq, the energy density behaves as p — > oo and therefore \p\ — > oo. Then the scale factor 
a is finite even if p ^ cx3. Therefore a > 1 case corresponds to type III singularity. 

• In a = 1 case, ii A > 0, there occurs the Big Rip or type I singularity but if A < 0, there does not appear future 
singularity. 

• In case 1/2 < a < 1, when t tQ, all of p, \p\, and a diverge if ^ > then this corresponds to type I singularity. 

• In case < a < 1/2, when t to, we find p, \p\ and a — > ag but 

Ina- \t-to\^^ . (100) 

Since the exponent (a — 1)/(Q! — 1/2) is not always an integer, even if a is finite, the higher derivatives of H 
diverge in general. Therefore this case corresponds to type IV singularity. 

• In case a < 0, when t to, we find p — > 0, a — > oq but \p\ oo. Therefore this case corresponds to type II 
singularity. 

Hence, the brief review of dark fluid FRW cosmology with specific homogeneous EOS as well as its late-time behaviour 
(singularities) is given (see jl^l for more detail). 

At the next step, we consider the inhomogencous EOS for dark fluid, so that the dependence from Hubble parameter 
is included in EOS. This new terms may origin from string/M-theory, braneworld or modified gra vity . One more 
motivation for such EOS comes from including of time-dependent bulk viscosity in dark fluid EOS |43 or from the 
modification of gravity. Hence, we suggest the following EOS 01 



P=-p + .f{p)+G{H) . (101) 



where G{H) is some function. 

As an example, one may consider the case 



f{p) = Ap^ -> !{p) + G{H) = -Ap^ - BH-"^ . (102) 

Here, it follows 

• In case a > \, for most values of /3, there occurs type HI singularity. In addition to the type HI singularity, 
when < /3 < 1/2, there occurs type IV singularity and when /3 < 0, there occurs type II singularity. 

• a = 1 case, if /3 > 1, the singularity becomes type III. If /3 < 1 and A > 0, there occurs the Big Rip or type 
I singularity. In addition to the type I singularity, we have type IV singularity when < /3 < 1/2 and type II 
when /3 < 1. 

• In case 1/2 < a < 1, one sees singularity of type III for /3 > 1, type Iforl/2</3<l (even for /3 — 1/2) or 
/3 = 1 and i?' > (i? > 0) case. In addition to type I, type IV case occurs for < /? < 1/2, and type II for 
/3 < 0. 

• In case a — 1/2, we have singularity of type III for /3 > 1, type I for 1/2 < /3 < 1 or /3 = 1 and B' > {B > 0), 
type IV for < /? < 1/2, and type II for /3 < 0. When /3==l/2or/? = 0, there does not appear any singularity. 

• In case 0<a<l/2, we find type IV for < /3 < 1/2, and type II for /3 < 0. In addition to type IV singularity, 
there occurs singularity of type HI for (3 > 1, type I for 1/2 < /3 < 1 or /3 = 1 and B' > {B > 0) case. 

• In case a < 0, there will always occur type II singularity. In addition to type II singularity, we have a singularity 
of type HI for /3 > 1, type I for 1/2 < /3 < 1 or /3 = 1 and B' > (B > 0) case. 

Thus, we demonstrated how the modification of EOS by Hubble dependent, inhomogencous term changes the 
structure of singularities in late-time dark energy universe. 

In general, EOS needs to be double- valued in order for the transition (crossing of phantom divide) to occur between 
the region w < —1 and the region w > —1. Then there could not be one-to-one correspondence between p and p. In 
such a case, instead of (|101|l . we may suggest the implicit, inhomogeneous equation of the state 

F{p,p,H)=0. (103) 
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The following example may be of interest: 



Here Cq and Hq are positive constants. Hence, the Hubble rate looks as 

16 



and 



P 



-P 



H 



4Hn ^ 



to) 



P 



(104) 



(105) 



(106) 



In H105(l . since t_ < <o < as long as <_ < t < <+, the Hubble rate H is positive. The Hubble rate H has a 
minimum H = Hq when t ^ to = (t_ + i+) /2 and diverges when t ^ t±. Then one may regard t ^ as a Big 
Bang singularity and t ^ t+ as a Big Rip one. As clear from (|106|l . the parameter w — p/ p is larger than —1 when 
t- < t < to and smaller than —1 when < t < t+. Therefore there occurs the crossing of phantom divide w — —1 
when t = to thanks to the effect of inhomogeneous term in EOS. 

In principle, the more general EOS may contain the derivatives of H, like H, H, ... More general EOS than H103|) 
may have the following form: 



f(p,p,h,h,h,---) =0 



Trivial example is that 



p^wp -H 



2^ 3(l + «;)^2 



(107) 



(108) 



By using the first and second FRW equations p = {3/k^)H^ and p = —{2/k^)H — {2,/k'^)H'^, Ea. (|108|l becomes an 
identity, which means that any FRW cosmology can be a solution if EOS H108f) is assumed. 



B. Origin of the inhomogeneous terms 

In order to give the support for an inhomogeneous EoS dark fluid |47j we will indicate how the extra terms may 
appear from the modifications of General Relativity. Let us start with the action of the scalar-tensor theory (|25|1 in 
sect^ The energy density p and the pressure p for the scalar field are given by 

p = ic^((^)02 ^ ^ p ^ icc>(0)02 _ v{c^) . (109) 

Since we can always redefine the scalar field as ^ P{4>) by an arbitrary function F{(p), we can choose the scalar 
field to be a time coordinate; (j) ~ t. Furthermore we consider the case that uj((j)) and V{(f>) are given by a single 
function /((/)) as (see sect lll D|) 

= -^f{<l>) , ^(0) = 4 + fW) ■ (110) 



Hence, one finds 

P= - 

Since p — f^^ {^Ky/pjij , we obtain the inhomogeneous equation of the state 



-/(0)', P = --/(0)' --/'(</-) ■ (111) 



p = -p--jir'i-xf^]] • (112) 



3, 

As an example, we consider the case 

/ = /o0' + /i , (113) 
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with constant /o and /i. Then (|112() gives 



/ 

If the matter contribution is neglected, from the first FRW equation, one finds 



= b + pf-%fl-4] • (114) 



t = ^P^H'. (115) 



Then we may rewrite (|114f) as 



0-(p + p)^-^J^(l-^) , (116) 



4/o fp (, fi 



k3 V 3 V H 

which has a structure similar to pU4|) . One may also consider the case that 



which gives 



(117) 



P=[-'^ + ^]p- (118) 



3/o 



Then we obtain dark fluid with homogeneous EoS with w = — l + 2/3/o. Note that despite the fact that specific scalar- 
tensor theory and dark fluid may predict the same FRW dynamics, their physical properties may be different. For 
instance, Newton law, number of FRW solutions, their stability properties may significally differ in two formulations. 

Another possibility to obtain the inhomogeneous EOS dark fluid is related with the modified gravity. As an 
illustrative example, the following action is considered: 

1 

2k2 



d^X^[—R + Crae.ttev + f{R)]. (119) 



In the FRW universe, the gravitational equations are: 

= -^H'^ + p- f (r^6H + 12H^ 

+6(h + - H-^j f (^R = 6H +12H'^^ , (120) 



(121) 



^ {2H + iH^^ +P + f {r = &H 

+2 (-H -iH^ + ^ + 2i74^ f (r = 6H+ 12H' 
\ dt^ at J \ 



Here p and p are the energy density and the pressure coming from £mattor- They may satisfy the equation of state 
like p = wp. One may now define the effective energy density p and p by 

(122) 



p = p- f (^R = 6H +12H^^ +6\^H + - H—j f (^R^6H + 12H'^^ 
p = p + f (r = 6H +12H'^ 



+2\^-H-3H^ + -^ + 2H—jf'{R = 6H+12H^j . (123) 

Thus, it follows 

p = wp+{l + w)f (r = 6H+ 12H^^ + 2 ^(-1 -3w)H -3{l + w)H^ + 

+ {2 + 3w)Hj\f'{R = m + \2H'^^. (124) 
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Especially if we consider the case that f = foR with a constant /o, we obtain 

p^wp + 6/oij + 6(1 + w)foH^ . (125) 

Thus, FRW dynamics (here, the same FRW equations) of some modified gravity may be equivalent to FRW dynamics 
in GR coupled with the inhomogeneous EoS dark fluid. Of course, the physics of two equivalent FRW descriptions 
(Newton law, perturbations structure, stability of FRW solution, number of FRW solutions, etc) may be quite different 
|2(TI |. In this respect, one should note that different modified gravities may also lead to the same FRW equations (see 
an explicit example in 13 )• Nevertheless, the details of FRW dynamics will be again different as modified gravity 
predictions should be compared with observations (for some attempts in this direction, see ,4£, ) . 

VII. DISCUSSION 

In summary, we presented the introduction to some classes of modified gravities considered as gravitational alterna- 
tive for dark energy. Specifically, the models which include f{R), f{G), f(R,G) terms or string-inspired GB-dilaton 
terms are considered (for recent review of higher derivative gravities, see [S^ and for discussion of their experimental 
manifestations, see |13). Special attention is paid to their versions where the action grows with decrease of the scalar 
curvature. In addition, the FRW universe filled with dark fluid with inhomogeneous EoS which may be induced by 
modified gravity is discussed. As it has been stressed in the introduction and as it has been explicitly demonstrated 
above, modified gravity is endowed with a quite rich cosmological structure: it may naturally lead to an effective 
cosmological constant, quintessence or phantom era in the late universe with a possible transition from deceleration 
to acceleration or (if necessary) the crossing of phantom divide. Moreover, some versions of modified gravity may 
pass the Solar System tests quite successfully. 

Nevertheless, we are still missing the modified gravity theory which, from one side, is as good as General Relativity 
in Solar System tests and from another side contains the fully consistent explanation for current cosmic acceleration. 
(Note that if universe is currently entering to phantom phase, then some deviations from General Relativity are quite 
expectable at large distances). It is quite possible that other versions of modified gravity: string- inspired gravity with 
fourth-order (or higher) in R corrections, theories which contain Traces of inverse Riemann tensor |5lj | or negative 
powers of squared Riemann and Ricci tensors 52] or even non-local Riemann and Ricci tensor terms (for instance, of 
the sort induced by Quantum Field Theory 53] with terms like R/ R^,yDR'^'^ , R/^i^R'^'^ / ROR, etc) or some combination 
of the above models with the ones under discussion should be considered as dark energy candidates as well. From 
another point, the next generation of precise observations should define the evolution of the cosmological parameters, 
and, perhaps, to rule out some of possible gravitational dark energy models leading to viable explanation of dark 
energy problem. 
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